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The formation of dibaryons with strangeness are discussed for the interior of neutron 
stars and for central relativistic heavy-ion collisions. We derive limits for the properties 
of H-dibaryons from pulsar data. Signals for the formation of possible bound states with 
hyperons at BNL's Relativistic Heavy-Ion Collider (RHIC) are investigated by studying 
j-C \ their weak decay patterns and production rates. 
O _ 

1. HYPERONS AND DIBARYONS IN NEUTRON STARS 

O 

There has been a lot of speculations about the appearance of new particles in the interior 
of neutron stars. A traditional neutron star consists of neutrons, protons, and leptons only. 
Cameron suggested first in 1959, that hyperons will also appear at high baryon density [|T|. 
A possible phase transition to quark matter was conjectured shortly after the quark model 
was introduced [01. Pion condensation M as well as kaon condensation were considered 
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later |4[]. Finally, strange stars, built out of absolutely stable strange quark matter, have 
been studied within the MIT bag model [f|||. The appearance of one or the other exotic 
phase in the core of neutron stars is still a matter of active debates. Nevertheless, recent 
developments in the field indicate, that hyperons are probably the first exotic particle to 
appear in neutron star matter at twice normal nuclear density. This result was consistently 
derived within various, different models: nonrelativistic potential models [0], quark-meson 
coupling model ||, relativistic mean- field models P-PL relativistic Hartree-Fock JT2[ , 



Brueckner-Hartree-Fock models 0,[TJj], and chiral effective Lagrangians |TS[. Hence, the 



onset of hyperons at 2po seems to be rather insensitive to the underlying details of the 
hyperon-nucleon interaction used. 

If there are a lot of hyperons in the interior of neutron stars, then it might be also 
possible to form dibaryons with strangeness. Besides the famous H-dibaryon proposed 



by Jaffe |TJ| built out of six quarks, the most recent version of the Nijmegen potential 
predicts the existence of deeply bound states of two hyperons, involving the S's and S's 
with maximum isospin ||17|| . Their results can be understood in terms of the underlying 



SU(3) symmetry for the baryon-baryon interactions (see P"T|). 

Let us assume in the following, that there exist bound states with hyperons. What 
would be the consequences for the properties of dense matter, as present in the interior of 
neutron stars? Bound systems in dense matter have been discussed before in connection 
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Figure 1. The equation of state of neutron star matter including a possible condensate 
of H-dibaryons. Uh denotes the optical potential of the H-dibaryon at normal nuclear 
density. For an attractive potential, the equation of state is considerably softened (taken 

from my 



with the liquid-gas phase transition of nuclear matter. They dissolve at rather low density 
due to effects in analogy to the Mott transition; deuterons disappear at l/15p and a 
particles around l/2p , as detailed e.g. in fL9 |. A similar transition will happen for 



bound states with hyperons. Hence, weakly bound systems will have no influence on the 
equation of state of neutron star matter — unless, there is a deeply bound state which 
can be treated as a quasiparticle. The H-dibaryon as well as heavier quark bags with 
strange quarks (strangelets) are such deeply bound states. They can actually be formed 
as precursor of the deconfmement phase transition in the core of neutron stars. The 
general condition for the appearance of a strangelet in s-wave with baryon number A and 
charge Z is that the effective energy of the particle equals its chemical potential 

E* s (k = 0) = M s + U(p) = Ax p n - Z x p e . (1) 

Here, U(p) stands for the effective potential felt by the strangelet at finite density. In 
order to estimate the effects of these exotic states, we focus on the lightest possible one, 
the H-dibaryon with A = 2 and Z = 0. We use an extended Relativistic Mean-Field 
Model which incorporates the H-particle in a consistent scheme. The Lagrangian part for 
the H-dibaryon field reads 

C H = D*^H*D^H - m* H 2 H*H (2) 

with a minimal coupling scheme for the scalar and vector meson fields 

Dp = + ig wH ■ (3) 
m* H = m H + 9aH ■ cr ■ (4) 
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For the coupling constant to the vector field, we choose the simple quark counting rule as 
guidance, i.e. we set g U H = 4/3^^. We fix the scalar coupling to an optical potential at 
normal nuclear matter density in the range 

U h (Po) = 9*h ■ <r{po) + g^H ■ V ( Po ) = -30 ... + 30 MeV . (5) 
Figure [1] shows the resulting equation of state when using the parameter set TM1 and 



including effects from hyperons (see |18|] for details). The equation of state turns out to 
be quite sensitive to the choice of the optical potential of the H-dibaryon. The phase 
transition to the H-dibaryon condensate is of second order in all cases studied. For an 
attractive optical potential, Uh(po) = —30 MeV, the equation of state has a plateau-like 
behaviour, i.e. the pressure is very slowly rising with energy density. 
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Figure 2. Neutron star constraints on the H-dibaryon from the mass of the Hulse- Taylor 
pulsar. A deeply bound H-dibaryon feeling an attractive potential in matter can be ruled 
out as the corresponding neutron stars have a too small maximum mass (from [fL8l). 



The drastic softening of the equation of state has crucial consequences for the maximum 
mass of a neutron star, which will be lowered below the acceptable value of the Hulse- 
Taylor pulsar of M = 1.44M . Hence, one can constrain the properties of the H-dibaryon 
in dense matter from pulsar data. These constraints are summarized in fig. [| as a function 
of the vacuum mass of the H-dibaryon and its optical potential at po- The calculations have 
been done by using two completely different parameterizations for the nucleon and the 
hyperon coupling constants (denoted as set TM1 and GM91 in the figure) and our findings 
seems to be rather insensitive to these choices. According to these calculations, the shaded 
area in the lower left side of the plot can be ruled out, as the corresponding neutron 
star equation of state gives too low a maximum mass. Therefore, a deeply attractive 
potential for the H-dibaryon with a mass close to the AA threshold is incompatible with 
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Table 1 

The hyperon weak decay amplitudes in SU(3) woa k compared to experimental data. All 



values are in units of 10 
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the measured Hulse- Taylor pulsar mass. Also, a H-dibaryon mass below the EN threshold 
requires even a repulsive potential in matter in our model calculation to get sufficiently 
heavy neutron stars. 

Finally, we remark that the existence of deeply bound states of two hyperons can have 
an indirect effect on the properties of neutron stars, though. If the hyperon-hyperon 
interaction is deeply attractive, it can give rise to a phase transition to hyperonic matter 
[ p0|1 - This phase transition can result in a drastically modified mass-radius relation for 
neutron stars, like neutron star twins with similar masses but smaller radii than ordinary 



neutron stars ||21|| . Hence, measurements of the mass and the radius of neutron stars will 
provide important insights into the equation of state and may reveal that there is indeed 
something strange going on at high density. 



2. STRANGE DIBARYONS IN RELATIVISTIC HEAVY-ION COLLISIONS 

There is the possibility to probe the hyperon-hyperon interaction in less remote places 
than neutron stars. Terrestrial relativistic heavy-ion collisions will produce dozens of 
hyperons in a single central collisions of two heavy nuclei, and for RHIC even close to 200 
hyperons are expected [p2| . The hyperons, sitting close in phase-space, can coalesce to 
form a bound state or a resonance which is decaying afterwards. The final products will 
be measured then in the detectors. 

The production rates for dibaryons with hyperons at RHIC has been estimated in the 
coalescence model to be in the range of 10~ 2 to 10~ 4 per single event. Effects from 
the details of the dibaryon wavefunction have been found to be rather small ||22|| . The 
rapidity distribution is rather flat, so that dibaryons are also produced at forward and 
backward rapidities where the decay length is considerably longer than at midrapidity. 
The production rates can be enhanced dramatically in two ways which are not taken into 
account in the above coalescence estimates. If a quark-gluon plasma is formed, strange 
quark- ant iquark pairs are produced more abundantly so that the total initial number of 
produced strange hadrons increases. If matter is created in a chirally restored phase, 
the effective hyperon masses are reduced drastically which will result in an enhanced 
hyperon-antihyperon yield. 
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Figure 3. Branching ratios for a bound AA state as a function of its binding energy. The 
nonmesonic decay mode AA — > X~+p dominates for a binding energy of 4 MeV or more 
(taken from fl22f). 



To detect strange dibaryons, their weak decay patterns have to be investigated. Start- 
ing point is the weak nonleptonic decay of hyperons in vacuum. Both, s-wave and p-wave 
amplitudes have been measured and correspond to a parity violating and a parity con- 
serving amplitudes, respectively. The standard approach for describing these amplitudes 
is by means of chiral perturbation theory. The p-wave amplitudes are beyond leading 
order and are usually derived within the pole model. The basic version fails to describe 
the p-wave amplitudes. This constitutes the so called s-wave/p-wave puzzle for the weak 
nonleptonic decay of hyperons (for a discussion and more elaborate ways to remedy the 
situation see p3] and references therein). We have found a simple way to parameterize 



both amplitudes, s-wave and p-wave, in terms of pure SU(3) symmetry for the weak in- 
teractions |[22|| . The effective Lagrangian involves the baryon octet B, the pseudoscalar 
octet P, and the Gell-Mann matrix X 6 , which ensures the AI = 1/2 rule and a change of 
hypercharge by one unit: 

C = D ■ TvBB [P, A 6 ] + F ■ TrB [P, A 6 ] B 

+ G ■ TrPP 75 PA 6 + H • TrPA 675 PP + J ■ TrP{P, A 6 } 75 P . (6) 

The first two terms give the s-wave contributions, while the other three the p-wave contri- 
butions. Table |I| compares the resulting amplitudes with the measured ones when setting 
D = 4.72, F = -1.62, G = 40.0, H = 47.8, and J = -7.1 (in units of 10" 10 ). We 
conclude that the above model as defined in eq. (^j) gives a reasonable description of all 
amplitudes. 

The weak decays of possible dibaryons with hyperons are then calculated by adopting 
a Hulthen-type wavefunction with varying binding energy. The nonmesonic decay is 
described by pion- and kaon- exchange diagrams, where the weak coupling constants are 
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Figure 4. Branching ratios for a bound S°p state versus its binding energy. The decay 
mode S°p^ Ap has the largest branching ratio for binding energies of 1.5 MeV or more 
(from H). 



given by eq. and the strong coupling constants by SU(6) symmetry. We note in passing, 
that the state-of-the-art calculations for the weak nonmesonic decay of hypernuclei rely 
on the (wrong) pole model (see p4[ and references therein). While this will not affect 



the pion exchange diagrams, accidentally, it will certainly alter the coupling constants for 
the kaon exchange diagrams. The use of the SU(3) symmetric model might be interesting 
to pursue in this direction to investigate possible effects, e.g. for the ratio of neutron- to 
proton-induced decay which comes out too small in the standard approach. 

Figure [3] depicts the branching ratio for a bound AA system versus the binding energy, 
and fig. f| shows the case for a possible S°p dibaryon. The mesonic decay channel domi- 
nates for binding energies of a few MeV. Then, for higher binding energies, the nonmesonic 
decay channels take over. For AA, the dominant nonmesonic decay mode is to a £™ and 
a proton, the same as for the hypothetical H-particle. The S°p decays mainly to a A and 
a proton for binding energies of just 1.5 MeV and more. This weak decay has the same 
decay topology as the weak decay of the S~ and the Q~ , and both particles have readily 
been measured in relativistic heavy-ion collisions p5| . Other interesting decay topologies 
occur for the doubly negatively charged dibaryons and all of which 



have been predicted to be bound in the recent Nijmegen model [17|. Here, a negatively 



charged object decays into two negatively charged tracks — a unique decay prong which 
should be easily seen in tracking devices. The calculated lifetimes of all dibaryons are 
lying just below the one for the free A with decay lengths in the range of cr = 1-5 cm. 

There are in principle three ways to detect possible short-lived strange dibaryons (for an 
experimental investigation we refer to the detailed feasibility study of Coffin and Kuhn for 
the case of the H-dibaryon at the STAR detector [^6|). Firstly, if the dibaryon is bound, 
one can look for exotic decays in a TPC, like a charged track splitting into two charged 
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ones. Secondly, if the mass of the strange dibaryon is close to threshold, it will show up 
in the invariant mass spectrum of its decay products by using background subtraction. 
Note, that this method is also sensitive to dibaryon resonances, if their decay widths are 
not too large ||27|| . Thirdly, the hyperon-hyperon interaction and possible broad dibaryon 
resonances can be probed in correlation functions as the low momentum part is sensitive 
to final state interactions 
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